Perpetual American options within CTRW's 
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Abstract 

Continuous-time random walks are a well suited tool for the description of market 
behaviour at the smallest scale: the tick-to-tick evolution. We will apply this kind 
of market model to the valuation of perpetual American options: derivatives with 
no maturity that can be exercised at any time. Our approach leads to option prices 
that fulfil financial formulas when canonical assumptions on the dynamics governing 
the process are made, but it is still suitable for more exotic market conditions. 
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Since their introduction in 1965 by Montroll and Weiss [I], continuous-time 
random walks (CTRW's) have been applied to virtually any field in which 



one wants to capture the smallest-scale properties of a given random system, 
and finance is not an exception. In fact, an intense activity in this area have 
been developed in recent years, as it can be inferred from the extensive (and 
exhaustive) list of references included in the review that Scalas [2J published 
in 2006. However, most of the work involving CTRW models was intended 
to reproduce empirical statistical properties, like probability density function 
(PDF) of returns [3], or the mean exit time of the process out of a given 
region [4f5] . The amount of such articles that deal with problems with some 
application on option pricing, as in Ref. [6], is not so large instead. 

In this article we present pricing expressions for several perpetual American 
options. The results were obtained by applying standard CTRW techniques 
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on a market model whose time evolution is step-like. The fact is that option 
pricing when underlying stock returns are discontinuous is not a new problem 
at all [7], and it is under active research nowadays in mathematical finance [8]. 
The main differences between these approaches are twofold: we do not assume 
either that the process evolution is diffusive between changes, CTRW collects 
all the stochastic (tick-by-tick) dynamics, nor transactions constitute neces- 
sarily a Poisson process. This might allow us to explore inefficient market 
models. When a risk-neutral market evolution can be defined we obtain fair 
option prices, from which classical Black-Scholes (BS) results can be recovered. 



2 The process 



We will model the market evolution of stock prices, S(t), through the logarith- 
mic return X(t) = \n[S(t)]. We will assume that the stochastic process follows 
the most common version of the CTRW: X(t) shows a series of random incre- 
ments or jumps at random times • • • , t_i, t ,t 1 ,t 2 , ■ ■ ■ , t n , ■ ■ ■ , the transaction 
times, remaining constant between these jumps. Therefore, the resulting tra- 
jectory consists of a series of steps as shown in Fig. [TJ We will assume that wait- 
ing times, At n = t n — t n _i, and random returns, AX n (t n ) = X(t n ) — X(£ n _i), 
are (mutually) independent and identically distributed random variables, de- 
scribed by two PDF's which we will denote by ip(t) and h{x) respectively. 

X(t) 
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Fig. 1. A sample trajectory of the X{t) process along with the corresponding value 
of the random variable r, the exercise time. 

A key magnitude in the CTRW formalism is the propagator, p(x — x , t — t ), 
the PDF of X(t) = x when one knows the value of the process right after a 
jump, X(t ) = x . Under our previous assumptions the propagator fulfils a 
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renewal equation [2] that can be readily solved in the Fourier-Laplace domain: 

P(U>,8) = - 1 

where tildes and hats stand for Fourier and Laplace transforms respectively. 



3 Option pricing 



Options are contracts between two parties, sold by one party to another, that 
give the buyer the right, but not the obligation, to buy (call) or sell (put) 
shares of the underlying stock at some prearranged price, the strike price K, 
within a certain period or on a specific date, the maturity or expiration time T. 
Therefore options are contingent claims with their present value determined 
by the discounted value of the expected profit under a risk-neutral measure: 

C{S ,t ) = E[P(S,r)e-^- to h T ^(t )], (2) 

where r is the risk-free interest rate, r is the actual exercise time, P(S, t) is 
the pay-off function, and T{t) condenses the available information up to time 
t. We will concentrate our attention on pay-off functions of the current value 
of the asset, like in the case of vanilla calls (+) and puts (-) where 

P ± (S,t)=±[S(t)-K]1 s{t ^ k , (3) 

and K is constant. Another broadly used pay-off is P ± (S, r) = 1s(t)>k i which 
corresponds to binary or digital call and put options. 

When the option can be exercised at the end of the contract lifetime solely, 
r = T, the option is said to be European. If the option can be exercised at 
any time before expiration it is called American, and r becomes a stochastic 
magnitude. Since the contract becomes worthless after maturity the option 
holder must find, under this constraint, the optimal exercise boundary, H(t): 
the stock price above (below) which it is better to exercise the call (put) than 
to keep the option alive. We will define the exercise time r ± as the first time 
the underlying crosses the threshold given that at present time, to, the spot 
price of the asset, So, lies in the proper side of the boundary: 

r± = min {t > t ; S(t) ^ H±(t)\S g H±} . 

Note that here the boundary function H(t) has to be assessed while solving 
the whole problem: the option price must always be greater or equal than the 
pay-off function, otherwise you can readily buy the option and execute it just 
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afterwards, thus obtaining risk- free earnings. Therefore -ff ± (t) must be settled 
in such a way that 

C ± (S,t)>P ± (S,t), (4) 
holds [9] . We must stress again that Eq. ([2]) leads to the right option price from 
the financial point of view if we use a risk-neutral measure. Under such market 
measure, and when the underlying stock pays no dividend, we must have 
E[y(5', t)e~ r< " t ~ t °'\!F{to)\ = Y(Sq, to), for any quoted quantity in the market. In 
particular ~E[C(S, t)e~ r ^~ to ^\!F(to)\ = C(So,to), which is true by construction, 
and ¥,[S(t)e~ r ( t ~ t °* ) \J-'(to)] = So- This martingale condition [8j for the risk- 
neutral evolution of the asset implies that our propagator must fulfil that: 

p{uj = —i, s + r) = s" 1 . 

If we put this expression into Eq. ([I]) we will arrive to the conclusion that 
waiting times ought to be exponentially distributed, ip(t) = \e~ xt , with 

r 

A = „ . (5) 

h(u = —i) — 1 

Let us analyse the consequences of our previous results. The physical measure 
and the risk-neutral measure may differ, but they must describe the same kind 
of process. This means that, if the number of actual market transactions does 
not follow a Poisson distribution, we will not be able to define a risk-neutral 
market measure, and the market will become inefficient, because then the 
CTRW has memory and it is not Markovian. Conversely, if waiting times have 
an exponential PDF, the risk-neutral market measure can be obtained after 
imposing Eq. (j5J). Note however that Eq. (J^J) is not a definition for A, but a 
constraint that involves all the physical parameters, as we will show below in 
a practical example. This is a typical feature of incomplete market models. 



4 Perpetual American options and survival probabilities 



Perpetual options, i.e. derivatives with T — > oo, are not actual traded con- 
tracts. However, they have practical interest because they represent the lim- 
iting value of a far-from-maturity contract, and therefore they may help in 
the pricing process if the theoretical price cannot be computed [10J — which 
is the most common situation for American options. Moreover, from a pure 
academic point of view, perpetual American options bring a major simplifi- 
cation to the issue: the optimal exercise boundary is constant H ± {t) = Hq, 
given that the problem is stationary. Then we have: 

C ± (So,to)=E[P ± (S;r ± )e- r ^ ± - t ^(U 

Observe that, since the process is discontinuous, the actual value of the stock 
price when we execute the option call (put) will be greater (lower) than Hq. 
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Therefore P ± {S; t^) is random in general. This is not the case for binary 
option prices, D ± (S ,t ), because P ± (S; r 1 * 1 ) = 1 and Hq = K : 

D ± (S ,t ) = E[e- r ^- to ^(t )} = 1 -r$| (a = r;x ) . 

Here we have introduced the complementary distribution function of r ± , 

®k (t ~ t ;x ) = Pr{t <t< r^Xito) =x = ln{S )} , 

which, within the CTRW formalism, is usually referred as the survival prob- 
ability (SP), S[ a ,b](t — to;xo). The SP measures the likelihood that a process 
beginning in a given region xq G [a, b] remains there all the time interval t — to- 
In all the cases we are dealing with, one of the boundaries of the region is in- 
finite. However, from a practical point of view, it is better to obtain the result 
for a finite region first, and compute the right limits afterwards. 

The Laplace transform of the SP fulfils by itself a renewal equation [6] that 
cannot be solved for a general kernel h(x). Therefore, we need to specify a 
functional form for the jump PDF, and we have focused our attention on the 
asymmetric two-sided exponential case, 



h(x) 



IP 
7 + P 



-px 



l x> o + e^ x l x 



x<0 



(6) 



This model is very suitable for our purposes for two main reasons. On the one 
hand, we can eventually recover the BS results in a certain limit. And, on the 
other hand, the solution to this problem is intricate but obtainable, 



S[a,b]{s;x ) = -<1 ^-e 



QX-() 



/3(xo-a-b) 



p ( 7 + 0) e -*._ 7 (p_0) e -/» 

7 (p - a)e aa - p( 7 + a)e ab 



where we have defined some auxiliary magnitudes: 



P-l . 
a = — h 



P + 7 



-7P^(s), /3 



P-l 



2 V V 2 J * i > \ / ~ 2 
A—(p — a)(j + (3)e^ )a - (p - /?)( 7 + a)e^ b . 



P + 7 



7PV'0), 



We must remember, however, that we only need the value of S^&^s; Xq) when 
either a — ► — oo or b — > +oo, what simplifies the final expressions: 



7 4>(s) /S 



s s 7 + a \Kq 



K ( S > X 0) 



1 p $( S ) /S \ p 
s s p — (3 \K ) 
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We can compare these results with the same outcomes for the Wiener process: 



♦W^) = | {i -(!)"}, *i(^) = ;{i-(|)'}. 

with 

a = + 4V^ 2 + 2a 2 s, = - -4v^ 2 + 2a 2 s. 

Here a is the volatility, the square root of the diffusion coefficient, and d is 
the drift of the diffusive process. The two formulas coincide for small values 
of the mean sojourn time p, for which we can expand ip{ s ) ~ 1 — ps, and for 
large values of p and 7, once one identifies ■& = (7 — p)/ r ypp, and cr 2 = 2/jpp. 



5 Risk-neutral prices 



In this section we will present results for the jump PDF in Eq. (J6]) when we 
have adopted a risk- neutral market measure, that is ip(t) = \e~ xt . We will 
proceed by replacing A according to Eq. (jSJ). This choice is feasible only when 
7 > p — 1 > 0, because we need h(u = —i) to be bounded, and A must be 
positive definite: 

7 - p+ 1 

The above conditions simplify the values of a = 1 and (3 = — (7 — p + 1), after 
we have set s — r, and lead to the following expressions for the live binary 
call (So ^ Kq) and put (Sq ^ K ) option prices|_j 

£ + (S ) = ^#, D-(So) ^- 1 ^o N7 " P+1 



P ^0 7 VS 

Note that in both cases the prices are discontinuous, since D + (So > K ) = 
D~(So < Kq) = 1. This feature disappears when considering continuous trad- 
ing, A — >• 00. We can approach to this limit by letting p — > 00 and 7 — > 00, 
but in such a way that the difference remains finite oo>7 — p + l = £:>0. In 
fact, in this case we recover the BS prices [HH2], since then a — a = 1, and 
/3 = ]f=-e= -2r/a 2 . 

We can now consider the usual vanilla pay-off which we introduced in Eq. ([3]). 
The problem here is that, as we advanced in Section HI -P ± (S; t^) is a random 
magnitude. When h(x) follows Eq. (jSj), however, K[P ± (S;T ± )\r :iz } does not 
depend on t ± . Then the vanilla option price l^ ± (So) can be easily computed: 



We will drop hereafter the dependence on to, because it is merely formal. 
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V + (S 
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V-(So) 
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(So > Ho), 



in terms of Hq . When we impose condition (j3J) on V + (Hq) we obtain an 
inconsistent expression for finite values of Hq and the option is never exercised, 
V + (S ) = S - When we demand V~(Hq) = (K — Hq), in turn, we get 



(7 + l)( 7 -p + l) 
7(7 ~P + 2) 



K <K. 



The continuous-trading limit leads again to the BS formula [mfT3] . as it can 
be observed in Fig. [2j There we represent put option prices for different values 
of p and 7. In order to reach the diffusive limit for a given choice of r and a, 
we have imposed the relationship 7 = p— l + 2r/a 2 . 
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Fig. 2. Option prices for different values of p. We show several plots of vanilla puts 
in terms of the moneyness Sq/K, after setting 7 = p — 1 + e, and e = 10. We have 
also represented the BS price for r = 5% and a = 10%. 



6 Conclusions 



We have argued for the convenience of the use of CTRWs in the modelling 
of stochastic processes in finance. CTRW is a well suited tool for representing 
market changes at the lower time scale, the transaction-to-transaction evo- 
lution. We have shown how CTRW-based results can be easily adapted to 
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financial terminology with little effort. In particular, we have analysed their 
possible applications in option pricing problems, by relating the computation 
of survival probabilities with the price of perpetual American derivatives for 
the two most ubiquitous flavours: binary options and vanilla options. 

We have considered risk-neutral scenarios for which we have obtained results 
that are fair from the financial point of view. Moreover, we have shown how 
classical BS results can be obtained from these expressions under certain lim- 
its. Finally, we have pointed out that our approach may be fruitful in the 
future in exploring a particular class of inefficient market models. There are 
multiple evidences [2"f4"] pointing to the fact that transactions are not exponen- 
tially distributed. In such a case the CTRW machinery may be very helpful 
in providing solutions for a general pausing time distribution. 
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